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Abstract 

We define a new average — termed the resolvent average — for positive semidefinite matrices. 
For positive definite matrices, the resolvent average enjoys self-duality and it interpolates be- 
tween the harmonic and the arithmetic averages, which it approaches when taking appropriate 
limits. We compare the resolvent average to the geometric mean. Some applications to matrix 
functions are also given. 
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1 Introduction 

Let Aj, i = 1, . . . , n be positive semidefinite matrices, \ > with Yl7=i ^« = 1 an d Id : — > M. N 
be the identity mapping. For 

A = (Ai,...,A n ), A = (Ai,...,A„), 

we define 

TZ^A, A) = [Ai(Ax + fi- 1 Id)- 1 + • • • + A n (A n + /i" 1 Id)- 1 ] ^ - pT 1 Id, (1) 
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and call it the resolvent average of A. This is motivated from the fact that when \x = 1 

(fti(A,A) + Id)~ 1 = Ai(A 1 + Id)~ 1 + --- + A n (A ri + Id)~ 1 , (2) 

which says that the resolvent of TZ\(A, A) is the (arithmetic) average of resolvents of the Ai, with 
weight A = (Ai, . . . , A n ). The resolvent average provides a novel averaging technique, and having 
the parameter [i in lZu(A,\) will allow us to take limits later on. We denote the well known 
harmonic average and arithmetic average by 

H(A,\) = (\ 1 A^ 1 + . . . + x n A^)-\ 

A(A, A) = Ai^i + • • • + X n A n , 

respectively. In the literature, (A± + ••• + A' 1 )' 1 is called the parallel sum of the matrices 
A 1 ,...,A n ; see, e.g., [D El EJ El Efl E21 [23] . 

The goal of this note is to study relationships among the resolvent average, the harmonic average 
and the arithmetic average of matrices. Our proofs are based on convex analytical techniques and 
on the proximal average, instead of the more commonly employed matrix diagonalizations. 

The plan of the paper is as follows. After proving some elementary properties of TZ^(A, A) in 
Section [21 we gather some basic properties of proximal averages and general convex functions in 
Section [3l The main results, which are given in Section HI state that 

H{A, A) * K^A, A) * A(A, A), 

lim Hu(A, A) = A(A, A), lim TZJA, A) = H(A, A), 

and that 7Z^(A, A) enjoys self-duality, namely [7Z^(A, A)] 1 = TZ^-i (A -1 , A). In Section[5l we show 
that the resolvent average and geometric mean have strikingly similar properties, even though they 
are different. 

Notation: Throughout, M> N is the standard iV-dimensional Euclidean space. For A > 0, 

J A = (Id+^)~ 1 , X A = \-\U-Jxa), (3) 

are called the resolvent of A and Yosida X-regularization of A. A function / : ~R N — > ]— oo, +oo] = 
R U {+00} is said to be convex if its domain is convex and 

/(Ax + (1-A)y) < A/(x) + (l-A)/(y) V x,y G R*,0 < A < 1, (4) 

with / being strictly convex if ((4]) becomes a strict inequality whenever x ^ y. The function / is 
proper if f(x) > — 00 Vx € M. N and /(xo) < +00 for some xq £ R . The class of proper lower 
semicontinuous convex functions from M> N — > ]— 00, +00] will be denoted by T. For / G T, df 
denotes its convex subdifferential: df(x) = {x* £ R N : f(y) > f{x) + (x*,y — x) My £ W N }. 
If / is differentiable at x, then df(x) = {V/(x)}. /* denotes its Fenchel conjugate given by 
(Vx* £ R*) f*(x*) = sup x {(x*,x) - f{x)}. For a > 0, a*f = af{-/a). If f,g £ T, fUg stands 
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for the infimal convolution of /, g given by (fdg)(x) = inf{/(xi) + 5(22) : x i + x 2 = %} Vx £ R^. 
When A : 1^ — ► R w is linear, the quadratic form q A : R N -► R is defined by 

ffil (a;) = ±(Aa;,a;) Vx G R w , 

and we also use gid = j interchangeably. For convex functions fi, ■ ■ ■ , f n , we write 

/ = (/l > • • • ) fn ) > / = (ft i • • • > f n ) • 

In the space of N x iV real symmetric matrices, (resp. §+ + ) denotes the set of N x N 
positive semidefinite matrices (resp. positive definite matrices). For X, Y € S N , we write Y <X 

if x - y g and y -< x if x - y e 



2 Basic properties 

In this section, we give some basic properties of 7£ M (A, A). 
Proposition 2.1 JVe /mwe 

Jf/R,fi(A,X) = H + ^nJnAn, (5) 

"(ft M (A,A)) = AiMi + --- + A n M n . (6) 

Proof. Multiplying ([1]) both sides by [i gives 

fj,n p (A, A) + Id = [Ai(M! + Id)- 1 + • • • + A n (M„ + Id)" 1 ] -1 . 
Then (0) follows by taking inverse both sides and using Q. 
By ([5]), we obtain that 

(Id - J^ M ( a, A) ) = A i ( Id - J Mi ) + ••• + A„ (Id - J Mn ) • 
Dividing both sides by fi, 

M~ X (Id -J^ A ,\)) = Aim -1 (Id - J Ml ) + • • • + A„/i _1 (Id -J^aJ- 
It remains to use ([3]). ■ 

Proposition 2.2 Let A = (Ai, A^ 1 , . . . , A m , A^ 1 ), A = (^, . . . , ^), and fx = 1. TTien 
ft„(A,A)=Id. 

Proo/. This follows from © and the identity (A + Id) -1 + (A' 1 + Id)" 1 = Id . ■ 
Proposition 2.3 Let A = (Ai, . . . , A\). Then K a (A, A) = A 1 . 



3 



Proof. We have 

U^A, A) = ((Ai + • • • + AnXAi + jT 1 Id)- 1 )" 1 - fj,' 1 Id 

= ({A 1 + fi- 1 Id)- 1 )' 1 - fild = A 1 + fi' 1 Id- fi- 1 Id = Ai, 
which proves the result. ■ 

Note that for A,P>£ we have 

A y B A~ l < B- 1 (7) 

and 

A y B ^ A' 1 ~< B~\ (8) 
see, e.g., [21 Corollary 7. 7.4. (a)] and [16j Section 16. E] or [101 page 55]. 

Proposition 2.4 Assume that (V i) A{,Bi E 8^ and >: Then 

n^(A,x)hn^B,\). (9) 

Furthermore, if additionally some Aj y Bj, then 1Z^(A, X) y 1Z^(B, A). 
Proof. Note that V fi > 0, 

+ /i" 1 Id ^ 5j + Id >- 0, 

so that 

■< (Ai + fi- 1 Id)" 1 ^ {B, + /i" 1 Id)- 1 , 
by ©. As 8+ and S^X + are convex cones, we obtain that 



o -< MA + m _1 id) -1 1 Yl Xi ( Bi + ^ Id ) _1 - ( 10 ) 

8=1 1=1 

Using ([7]) on (jTUJ) , followed by subtracting /i- 1 Id, gives 

n n 

[ h{Ai + pT 1 Id)- 1 ] " X - iT 1 Id b [ HBi + pT x Id)" 1 ] _1 - ^ Id, 



- 1 ] " 1 — M 1 Id >r [E x ; ° ' ■■- iT - ,1 - 11_ ' 1 
i=i i=i 

which establishes Q. The "Furthermore" part follows analogously using 



Theorem 2.5 Assume that (V i) Ai G SV. Then Ttn^A, A) € §+. Furthermore, if additionally 
some Aj G S^ + , i/ien TZ^A, A) G S^ + . 

Proof. This follows from Proposition 12.41 (with each l?j = 0) and Proposition 12.31 ■ 
We end this section with a recursion formula that may be verified directly using the definitions. 
Proposition 2.6 (recursion) We have 

TZn(A\ i ■ ■ ■ ! A n ; Ai , . . . , A ra ) = 7Z^ ( IZ^ (A\ , . . . , A n —i ; ^_\ n , • • • ; j ) , A n ; 1 — A n , A n 



3 Auxiliary results and facts 



The key tool in this note is the proximal average of convex functions, which finds its roots in 
[H UHl [20] , and which has been further systematically studied in (HJ [U [9] . 

Definition 3.1 (proximal average) Let (Vi) /, G T. The X-weighted proximal average of f = 
(/l) • • • 3 fn) with parameter fi is defined by 



p M (/, A) = f Ai(/i + ij)* + A 2 (/ 2 + Ji)* + • • • + An(/n + Jj)*) " Jj- 



(11) 



The function A) is a proper lower semicontinuous convex function on IR^, and it inherits 

many desirable properties from each underlying function ff, see [HE]. A fundamental property of 
proximal average is: 

Fact 3.2 (0 Theorem 5.1]) (p^(f, A))* = p^-i (/*, A). 

To give new proofs of Fact 13.41 and Fact 13.51 below, we shall need reformulations of Pu(f, A). 

Proposition 3.3 Let fi,...,f n ET and Ax, . . . , X n > with Ym=i Ai = 1. Then for every x € M n ; 
Pn(f, A)(x) 

= inf {Ai(/i + ii)(fi) + --- + A n (/ n + -i)(^))-ij(^) (12) 

xiH hr„=x fl Ai fj, X n J /i 



inf < 

xiH hx n =x 



1 n n 

Ai/i(fi) + • • • + A n/ „(^) + ^ £ £ A ^n? - ?n 2 ( 13 ) 



Ai A n 4/i t - ^ f Aj A,- 



inf I Ai/i(yi) H h X n f n (y n ) + ~[Axj(2/i) H H X n j{y n ) - j(Xiy x H h A n y n )] 

AiJ/H hA n j/,i=a; /i 

(14) 



inf 

AiyH hA 



{ 1 n n 1 
Ai/iG/i) + ■ • • + A n / n (y n ) + — EE AiA jll^-^'H 2 r ( 15 ) 
^ i=i i=i J 

inf (Ai/x^) + • • • + An/ n (^) + -[Axi(x - ^i)) + . . . + X n j(x (16) 

••+x n =x |^ Ai X n H Ax A n J 



xi H hx n =a; 

Furthermore, the infimal convolutions in il2\) -[T5 \) are exact 
Proof. Indeed, as 

(fi + 



fi + -j) =/;□(«), 
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it is finite-valued everywhere, we write 

/ = Ai * (A + -j)D ■ ■ ■ DX n ★ (f n + -j) - -j, 

by |25} Theorem 16.4]. That is, for every x, 

f(x) = inf (a^ + + ■■■ + X n {f n + -j)(^) : *i + • • • + x n = x) - -j(x), 

and the infimum is attained. Hence (|12p holds. 
Now rewrite (1121) as 



inf (ax/x^) + • • • + A„/ n (^) + -[Axi(^) + • • • + A n j(^) - j(x! + • • • + x n )] 

(17) 

= inf |Ai/i(yi)H V X n f n (y n ) + — [Xij(yi) H h X n j(y n ) - j{X\yi H h A n y n )] 

Aij/iH hA„j/ n =a; |^ /i 



Thus, p3j) - p5|) follow by using the identity 



Observe that 



i=l i=l i=l j = l 



Aii(xi H h x n - -i) H h A n j(xi H \-x n - 

= Aii(^-) H h A n j'( y 1 ) - H h x n ), 

Ai A n 



we have (fT6l) by (fT7]) . 

Fact 3.4 (0 Theorem 5.4]) (Ai/f + • • • + A n /*)* < p M (/, A) < Ai/i + • • • + X n f n . 



Proof. This follows from (|13p or (|15p . ■ 
Fact 3.5 ([3 Theorem 8.5]) Lei x € T/ien t/ie function 

]0,+oo[ — > ] — oo,+oo] : n i— > Pn(f ', A)(x) is decreasing. (18) 
Consequently, lim^ + p^f, A)(x) anc? lim M ^ +00 p^(/, A)(x) exist. In fact, 

lim p fl (f,X)(x) = supp^(/, A)(x) = (Ai/i H h A n / n )(x) (19) 



and 



lim ^(/,A)(x) = inf p M (/,A)(x) = (Ai */!□•• •□A n */ n )(x). (20) 



Proof. flUD follows from (PJ. $2$ also follows from (JTSJ). 
To see CED, by CE]), Vx G R^, 

pJf,X)(x) > Aiinf ( /i(xi/Ai) + -j'(x - ari/Ai) J H h A„inf ( f n (x n /K) + - x n /X n ) ) 

Xl \ / £n \ /X J 

= Aie^/i(x) H h A n e M / n (x), 

where e^/i = fiD(l/fj,j). Then 

Aie M /i H h Ane^/n < p^(/, A) < Ai/i H h A n / n , 

so that 

lim A) = Ai/i H h A n / n , 

since lim^ + e^fi = fi by |26|, Theorem 2.26 and Theorem 1.25]. ■ 

Fact 3.6 ( \25\ Theorem 25.7]) Let C be a nonempty open convex subset of R"^ , and let f be a 
convex function which is finite and differ entiable on C. Let /i,/2, ••• , be a sequence of convex 
functions finite and differ entiable on C such that linij^oo /j(x) = f(x) for every x E C. Then 

limV/i(x) = V/(x), VxG<7. 

i 

Ln fact, the sequence of gradients V/j converges to V/ uniformly on every compact subset of C. 
Fact 3.7 ([251 page 108]) Lei Q G §^ + . T/ien (g Q )* = q Q -i. 

Fact 3.8 ([251 Theorem 23.5]) Let f : R^ — > ]-oo,+oo] 6e a proper lower semicontinuous convex 
function. Then df* = {df)~ l . 



4 Main results 

We start by computing the proximal average of general linear-quadratic functions thereby extending 
[71 Example 4.5] and P Example 7.4]. 

Lemma 4.1 Let A{ G 6j G R^, rj G R. If each fi = qA t + (bi, ■) + r^, i.e., linear- quadratic, 
then Vx* ; 

n n 

p M (/, A)(x*) = 47MAA)(z*) + (x*, ( £ \i{Ai + /i- 1 Id)- 1 ) 1 + ^ Id)- 1 6 4 )+ 

i=l i=l 

n n 

(21) 
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In particular, if (Vi) fi is quadratic, i.e., bi = 0, r« = 0, then p^,(f, A) is quadratic with 

If (Vi) /i is ajffme, i.e., ^ = 0, then p^f , X) is affine. 

Proof. We have /j + = g^+^-i id) + {h, ') + r i and by Fact [377] 

(/i + M -1 j)*0O = <HAi+fi-i Id)-i(%* - bi) - n 

= Q(Ai+n-i id)-i - ( x *, ( A i + M" 1 Id) _1 &i) + 9(^+^-1 M )-i - n- 

Then (Ai(/i + M _1 j)* + • ■ ■ + V/* + M -1 i)*)0O = 

^AiU (Ai+At -i Id) -i(x*) - + ) u - " 1 Id) _1 6 i ) + 9(^+^-1 id)-i(^) -n) 

i=i ^ ' 

n n 

= ^E- l =i Ai^+^w)- 1 ^*) - (x*,J2 X i( A i + A* -1 Id) -1 &i) + ^ Ai(g {Ai+M -i Id) -i(6i) - ri). 

i=l i=l 

It follows that p fi (f,X)(x*) = 

n n 

^ELiMA+M- 1 ^)- 1 ]- 1 ^* + 5^A i (i4i + /i~ 1 Id)~ 1 6 i ) - ^ Ai (^ r (A l +^- 1 H)- 1 (^) - r - id ^ 

i=l i=l 

As 

n 

?E£=i ^(Ai+M-i id)- 1 ]-! ( x * + X] A ^' + Id ) _1 ^) 

i=i 

n n 

= 3E2=! ^(^+^-1 ld)-i]-i (a?*) + (a?*. E A »^ + ^ Id ) _1 ] -1 E A * ^ + ^ Id ) _1 ^) + 

1=1 8=1 

n 

?E2=i M^+M" 1 Id)- 1 ]- 1 !^ A ^ + ^ Id ) _1 ^)' 

8=1 

we obtain that p^if, A)(x*) = 

n n 

Q\ELi MAi+n- 1 h)-i]-i- m -ih(^) + (x*, E Ai(A + a^ 1 Id)" 1 ]- 1 E Ai(^ + /x- 1 Id)-\)+ 



8=1 1=1 



?EteiA ( (>li+A*-lId)-i]-i(J^^t(>*« + ^ ljd ) 1& - E Ai (^+M- 1 Id)-i(&i) 

1=1 8=1 

which is (|21[) . The remaining claims are immediate from (|2ip and that 1Z^(A,\) = when 
(V i) Ai = by Proposition [273J ■ 

We are ready for our main result: 



Theorem 4.2 (harmonic-resolvent-arithmetic average inequality and limits) 

Let At,..., An € §^ + . We have 

(i) 

H(A, A) * n^A, A) * A(A, A); (22) 
In particular, TZ^A, A) G §>++• 

(ii) 7*^(A, A) -> A(A, A) uAen /j -> 0+. 

(iii) 1Z^(A, A) — > H(A, A) u>/ien /i — > +oo. 



Proof, (i) . According to Fact [37 

(AiA* + • • • + A n /:)* < Pm(/> A) < Ai/i + ■ • • + A n / n . (23) 
Let /j = q^. Using ((7aJ* = (by Fact I3.T|) and Lemma ETTI we have 



i-i , 



(AlJT + • • • + XnfnT = (Al^ 1 + • ' ■ + An^-0* = (^-1+...+;^- 

= 9(A 1 Ar 1 +...+A n ^ 1 )-l = 8W(AA)- ( 24 ) 

Ai/i H h A n / n = q\ 1 A 1 +-+x n A n = qA(A,\), (25) 

Pn(f,X) = QR lt (A,\y ( 26 ) 

Then (|23|) becomes 

Qh(a,\) < Qn„(A,x) < Qa(A,x) ■ 

As q x < qy X < Y, ((22j) is established. Since Aj G S^ + , A' 1 G §^ + , AiA^+- • • + A n ^" 1 G 8^ + , 
we have H(A, A) = (A1A7 1 + • • • + A^ 1 )- 1 G §^ + , thus TZ^A, A) G by ||22J). (Alternatively, 
apply Theorem 12.51 ) 

(ii) and (iii): Observe that (Vi) (\i*fi)* = \-if* = Xiq.-i has full domain, by |25^ Theorem 16.4], 

i 

(Ai/r + • • • + \ n f:r = (x x */!□••• da„ * /„) . 

By Fact [331 Vx G one has 

lim A)(x) = (Ai/i H h A n / n )(x), 

ft-t-0+ 

lim p M (/, A)(x) = (Aa/{ + • • • + X n f*)*{x). 

fl— > + oo 

Since (Vi) /j, /* are differentiable on M. N , so is p^if, A) by Corollary 7.7]. According to Fact 13.61 
Vx 

lim Vp M (/, A)(x) = AiV/i(x) + • • • + A„V/„(x), (27) 
/u— >0+ 

lim V^(/,A)(x) = V(A 1 /r + --- + A n /:)*(x). (28) 
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Moreover, the convergences in (|27p - (|28p are uniform on every closed bounded subset of M. N . Now 
it follows from that Vp M (/, A) = TZ^A, A), V(Ai/i + • • • + X n fn) = A(A, A), V(A 1 / 1 * + 

. . . + \ n f*)* = H(A,X). dSZD-dSHD transpire to 

lim TZ^A, X)x = A(A, X)x, (29) 

lim KJA,X)x = H(A,X)x, (30) 

where the convergences are uniform on every closed bounded subset of M . Hence (ii) and (hi) 
follow from CUD and ([30]). ■ 

Note that in Theorem I4.2f ii).(iii). there is no ambiguity since all norms in finite dimensional 
spaces are equivalent. 

Definition 4.3 A function g: D — > S , where D is a convex subset of §> N , is matrix convex if 
VAx,A 2 e D, VA g [0,1], 

g(XA 1 + (1 - X)A 2 ) ± XgiAx) + (1 - A) 5 (A 2 ). 
Matrix concave functions are defined similarly. 

It is easy to see that a symmetric matrix valued function g is matrix concave (resp. convex) if 
and only if Vx £ IR^ the function A i— > <7 g M) (^) is concave (resp. convex). Some immediate 
consequences of Theorem 14.21 on matrix- valued functions are: 

Corollary 4.4 Assume that (Vi) £ and 2^=1 ^« = ^ ™^ ^* ^ 0- 

(AiAi + • • • + A^)' 1 ^ A^ 1 + • • • + A n A~ x . 

Consequently, the matrix function X *— > X^ 1 is matrix convex on §+ + . 

Proof. Apply Theorem I4T21 equation ([22|) for A = (A^ 1 , ■■■ , A^ 1 ). ■ 
Corollary 4.5 For every \i > 0, i/ie resolvent average matrix function A i— > IZ^A, A) given by 

{Ax, ■ ■ ■ , An) ^ [Ai(Ai + /x- 1 Id)- 1 + • • • + A n (^ n + Id)- 1 ]- 1 - /i -1 Id 
is matrix concave on x • • • x §+ + . (31) 

For eac/i A = (Ai, • • • , A n ) wi£/i Aj = 1 and Aj > Vi, i/ie harmonic average matrix function 
(A\, • • • , ^4 n ) i — > (Ai-A^ 1 + • • • + \ n A~ )~ is matrix concave (32) 
on §+ + x • • • x . Consequently, the harmonic average function 

(xi, ■ ■ ■ ,x n ) i— > — — is concave (33) 

+ ■ • • T" X n 

on K ++ x • • • x 



10 



Proof. Set fi = qAi- Then Vx € M. N , we have p^(f, A)(x) = 

min ( (Xiq Al (xi)-\ \- X n qA n (x n )) + ^iqid(xi) -\ h^ _1 A n gi d (x n )) ) -//"Vdfa 

Ai^iH — hA„x„=x y y 

Since for each fixed (xi, . . . , x n ), 

(Ax, ■■■ ,A n ) i-> (AigAi^i) H h KqA n (x n )) + (^ _1 gid(zi) H h M _1 <?id(zn)) , 

is affine, being the infimum of affine functions we have that Vx the function 

(Ax,-- - ,A n ) ^ p^(f,X)(x), 

is concave. As p fl (f, X)(x) = qn IJ ,{A,\)(x) by Lemma l4~Tj this shows that Vx £ f w the function 

A = (Ax, ■■■ ,A n ) ^ qn„(A,\)(x) is concave, 

so Ah TZn(A, A) is matrix concave. 

Now by Theorem I4.2f iii). TZ^(A, A) — > 7i(A, A) when /x — > +co. This and (f3"Tj) implies that 

is also matrix concave. (j33|) follows from (|32p by setting iV = 1 and Ai = • • • = A n = 1/n. I 



Remark 4.6 Corollary 14.41 is well-known, cf. |26^ Proposition 2.56 on page 73]. Corollary 14.51 (|33p 
is also well-known, cf. Exercise 3.17 on page 116]. 

We proceed to show that resolvent averages of matrices enjoy self-duality. 

Theorem 4.7 (self-duality) Let (V i) A, L G and fj, > 0. A ssume that Yl?=i Ai = 1 with 

Ai > 0. Then 

[R-n(A, A)] 1 = IZ^-i (A -1 , A), i.e., (34) 



Ai(Ai + iT 1 Id)" 1 + • • • + A n (A„ + ^ U)- 1 ^ - fi' 1 Id 
Ai(Ai 1 + /xld)- 1 + • • • + A n (A" 1 + fild)- 1 ^ 
In particular, for \x = 1, [1Zx(A, A)] -1 = lZx(A~ l , A). 



/j Id. 



Proof. Let = g^. By Fact I3.2| (^(/, A))* = p„-i(/*,A), taking subgradients both sides, 
followed by using Fact 13.81 we obtain that 



d(p,(f,\)Y = (dp^f, A)) -1 = d(p^-i(f*, A)). 
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By Lemma El p^(/,A) = g^(A,A), V 1 A ) = ^ ^(A- 1 ,*)* we have 

dp lt (f,\)=K lt (A,\), 
dp tl -i(f*,X)=K fl - 1 (A-\\). 

Hence 

[T^A)]- 1 ^-!^- 1 ^), 
as claimed. ■ 

Remark 4.8 Although the harmonic and arithmetic average lack self-duality, they are dual to 
each other: 

[H(A, A)]- 1 = X X A^ + ■■■ + XnA- 1 = A(A-\\), 
[A(A, A)]- 1 = [MAj; 1 )- 1 + • • • + XniA' 1 )- 1 ] 1 = U{A-\\). 



5 A comparison to weighted geometric means 

If A, B £ the geometric mean is defined by 

A\B = AV^A-WBA-Wf^A 1 ' 2 . 

In general, the geometric mean of A±,...,A n £ for n > 3 is defined either as the limit of 
an inductive procedure or by the Riemannian distance without a closed form [21 [24"1 [19], [T5] . To 
compare the resolvent average with the well-known geometric mean, we restrict our attention to 
non- negative real numbers (lxl matrices). When A = x = (x\, . . . , x n ) with Xi € R+ and fi = 1, 
we write 

K(x, A) = Tl„(A, A) = (Ai(x! + l)" 1 + • • • + X n {x n + l)" 1 ) -1 - 1, 
and x~ 1 = (1/xi, . . . , l/x n ) when (Vi) Xi € 

Proposition 5.1 Let (V i) xi > 0, > 0. W^e /iawe 

(fy (harmonic-resolvent-arithmetic mean inequality); 

(Aix^ 1 H + XnX' 1 ) _1 < K(x, A) < Ai^i + • • • + X n x n . (35) 

Moreover, 1Z(x, A) = AiXi + • • • + A n x n ?/ and only if x% = ■ ■ ■ = x n . 
(ii) (self-duality); [R(x, A)]" 1 = K{x-\ A), 
(^mj If x = (xi, . . . , xi), i/ien 72- (a;, A) = x\. 

(iv) Ifx = (xi,x^ 1 ,x 2 ,X2 1 , . . .,x n ,x~ x ) and A = . . . , i), t/ien TZ(x,X) = 1. 
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(v) The function x \— > 7l(x,\) is concave on R++ x • • • x 

(vi) If x y y, then TZ(x, A) > TZ{y, A). 

Proof, (i): For ([35]) . apply Theorem I4.2f i) with /j, = 1. Now 7£(a:,A) = A1X1 + ••• + A n x„ is 
equivalent to 

(Ai(xi + l)- 1 + • • • + X n {x n + l)" 1 )" 1 = Aixi + • • • + \ n x n + 1, (36) 

As Yli=l Ai = 1, (f36|) is the same as 

1 x 1 1 

A 17 ^TT + --- + ^n- 



(xi + 1) (x n + l) Ai(xi + 1) + --- + A n (x n + 1)' 

Since the function x i— ► 1/x is strictly convex on we must have xi = • • • = x n . 

(ii): Theorem 14.71 (iii): Proposition 12.31 (iv): Proposition I2.2L (v): Corollary 14.51 (vi): 
Proposition 12.41 ■ 



Recall the weighted geometric mean: 

G(x,X) = x^x^---x^. 
Q{x,X) always has the following properties: 
Fact 5.2 Let (V i) Xi > 0, yi > 0. We have 

(i) (harmonic-geometric-arithmetic mean inequality): 

(XixJ 1 H h AnX" 1 ) 1 < Q(x, A) < Aixi H h A n x n . 

Moreover, Q(x, A) = Aixi + • • • + A n x n «/ and on/?/ cci = • • • = x n . 
(ii) (self-duality): [Q(x, A)] -1 = Q(x^ 1 , A), 
(mj If x = (xi, . . . , xi), i/ien C/(a;, A) = x\. 

(iv) If x = (xi,x^ x ,x 2 ,X2 , ■ ■ ■ jXnjX" 1 ) and A = (i, . . . , i/ien £7(x, A) = 1. 
fuj T/ie function x i— > £7 (a;, A) is concave on M ++ x • • • x 
(mj Ifx>zy, then Q(x, A) > £7(y, A). 

Proof, (i): See [25l page 29]. (ii)-(iv) and (vi) are simple, (v): See [261 Example 2.53]. 
The means 7£(a:, A) and (?(a:, A) have strikingly similar properties. Are they the same? 
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Example 5.3 (i). Let A = (|, |). When x = (0,1), Q(x,X) = but K{x,X) = ± so £(a;,A) / 
Tl(x,X). 

(ii). /s rig/ii £/ia£ ^(a;,A) < 1Z(x,\) V x G R+ + ? The answer is also no. Assume to the 
contrary that Q(x,X) < 7Z(x,\), V x € K++ x Taking inverse both sides, followed by 

applying the self-duality of G(x, A), lZ(x, A), gives 

^(a;, A)" 1 > K(x, A)" 1 = K( X -\ A) > ^(a;- 1 , A) = G(x, A) -1 , 

and this gives that ^(x, A) -1 = ft (a;, A) 1 so that G(x, A) = ft(a;, A). This is a contradiction to 
(i)- 

Finally, we note that the resolvent average can be defined for general monotone operators and 
that Theorem 14.71 holds even when A\, . . . , A n are monotone operators (not necessarily positive 
semi-definite matrices), in that situation one needs to use set-valued inverses. This and further 
details on the resolvent average for general monotone operators will appear in the forthcoming 
paper [5]. 
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